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"^H ! We re-examine the generation of the baryon asymmetry in the minimal supersymmetric 

standard model (MSSM) during the electroweak phase transition. We find that the dominant 
source for baryogenesis arises from the chargino sector. The CP-violation comes from the 
complex phase in the /i parameter, which provides CP-odd contributions to the particle 
dispersion relations. This leads to different accelerations for particles and antiparticles in 
the wall region which, combined with diffusion, leads to the separation of Higgsinos and 
their antiparticles in the front of the wall. These asymmetries get transported to produce 
perturbations in the left-handed chiral quarks, which then drive sphaleron interactions to 
create the baryon asymmetry. We present a complete derivation of the semiclassical WKB 
formalism, including the chargino dispersion relations and a self-consistent derivation of the 
diffusion equations starting from semiclassical Boltzmann equations for WKB-excitations. 
We stress the advantages of treating the transport equations in terms of the manifestly gauge 
invariant physical energy and kinetic momentum, rather than in the gauge variant canonical 
variables used in previous treatments. We show that a large enough baryon asymmetry can 
be created for the phase of the complex /^-parameter as small as ~ 10"^, which is consistent 
with bounds from the neutron electric dipole moment. 



1 Introduction 



It is a fascinating possibility that the baryon asymmetry of the universe (BAU) may have 
been generated at the electroweak epoch (for reviews, see The great attraction of this 
idea is that, in contrast to other mechanisms operating at higher energy scales, it involves 
physics which is being searched for at accelerators now. An a priori calculation of the baryon 
asymmetry, as accurate as that of the abundance of the light elements in nucleosynthesis, 
may still be unattainable at present, but we should nevertheless strive to compute it as 
carefully as possible. One hopes thereby to reach a definitive conclusion as to the feasibility, 
at least, of generating the BAU at the electroweak scale. 

While there are many theoretical motivations for considering extensions of the standard 
model (SM), in the present context we are also prompted to do so for the simple reason that 
the SM by itself appears unable to produce the observed BAU. The smallness of the CP vi- 
olation in the KM matrix appears to be in itself an insurmountable obstacle to baryogenesis 
in the SM (although there has been considerable debate on this subject [|]), and has moti- 
vated many studies of baryogenesis in extended models with additional CP violation leading 
to more efficient baryon production. In addition to this problem, moreover, the SM fails 
badly with respect to the sphaleron wash-out boundQ. Lattice studies have shown that for 
any value of the higgs mass, even well below the present experimental lower limit, the phase 
transition would be so weak that sphaleron interactions remain in equilibrium in the broken 
phase of the electroweak sector, causing the baryon asymmetry to relax back to essentially 
zero immediately after its generation [Q. 

Several extensions of the SM have been considered to overcome the sphaleron wash-out 
bound by strengthening the phase transition ||^, ^, 0, ^. Best motivated from the particle 
physics point of view is the minimal supersymmetric standard model (MSSM). Several recent 
perturbative and nonperturbative studies of the properties of the phase transition in this 
model H, |10|, O, [H] have shown that in a restricted part of the parameter space, the 
sphaleron bound can be satisfied. An important question is therefore whether for these same 
parameter values the generation of the observed BAU is possible. 

Baryogenesis in the MSSM has already been studied in several papers [|^, |T5|, |16 



^As discussed in this bound is predicated on the assumption that the Universe is radiation dominated 
at the electroweak epoch, and can be significantly weakened in non-standard (e.g. scalar field dominated) 
cosmologies. 



1 



T7| , [T8|, [T^, ^ ^ . The overall framework of the baryogenesis mechanism is essentially 



agreed upon: bubbles nucleate at a first order phase transition and the expanding bubble 
walls propagate through the hot plasma, perturbing the quasiparticle distributions from 
equilibrium in a CP-violating manner. Incorporating the effects of transport leads to a 
local excess or deficit of left-handed fermions over their antiparticles on and around the 
propagating bubble walls. This drives the anomalous baryon number violating processes to 
produce a net baryon asymmetry, which is swept behind the bubble wall where it is frozen 
in (assuming that the sphaleron bound is satisfied). Moveover, common to all methods is 
reducing the problem to a set of diffusion equations coupling the sourced species to the 
species that bias the sphalerons. These are coupled equations which have the general form 

AC + v^^' + r,(e. + ^j + ---) = s,, (1) 

where i labels the particle species and ^ = /ij/T is its chemical potential divided by temper- 
ature. Primes denote spatial derivatives in the direction (z) perpendicular to the wall, 
is the wall velocity, Fj is the rate of an interaction that converts species i into other kinds 
of particles, and Si is the source term associated with the current generated at the bubble 
wall. There is little controversy about the form of these equations, but little agreement exists 
as to how to properly derive the source terms S^. There are many different formalisms for 
obtaining the sources |^, but so far little effort has been made to see how far they 
agree or disagree with each other. We shall comment on this issue briefly in our conclusions. 

Here we shall use the 'classical force' mechanism (CFM) for baryogenesis 0, p8| , pO| , 
PT| . The CFM makes use of the intuitively simple picture of particles being transported 
in the plasma under the influence of the classical force exerted on them by the spatially 
varying Higgs field condensate. We assume that the plasma in this bubble wall region can 
be described by a collection of semiclassical quasiparticle states which we shall refer to as 
WKB states, because their equation of motion is derived using the WKB approximation 
expanding in derivatives of the background field. The force acting on the particles can be 
deduced from the WKB dispersion relations and their corresponding canonical equations of 
motion. This is a reasonable assumption when the de Broglie wavelength of the states is 
much shorter than the scale of variation of the bubble wall, i.e. X <^ i^, which is satisfied 
in electroweak baryogenesis; in the MSSM, the wall widths are typically ~ 6 — 14/T 
1^ , p6|] , whereas for a typical excitation A ~ 1/T. Given these conditions one can write a 
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semiclassical Boltzmann equation for the distribution functions of the local WKB-states 



{^t + ^rg.^^ + F■ 9p)/, = C[f„ f„ ...]. (2) 

where the group velocity and classical force are given respectively by 

Vg = dp^u] F = p = uJYg. (3) 

Here Pc is the canonical, and p = uj\g the physical, kinetic momentum along the WKB 
worldline. Note that we treat the transport problem here in the kinetic variables - in which 
the Boltzmann equation has the non-canonical form of (|^) - rather than in the canonical 
variables used in previous treatments. As will be discussed in more detail below, this choice 
has the simple advantage of circumventing all the difficulties associated with the variance 
of the canonical variables under local phase ('gauge') transformations of the fields in the 
Lagrangian. In these kinetic variables it is also more manifestly (and gauge independently) 
clear how, because of CP-violating effects, particles and antipartices experience different 
forces in the wall region, which leads to the separation of chiral currents. The explicit form 
of Vg and F in a given model can be found from the WKB dispersion relations, as we will 
illustrate in sections 2 and 3. The Boltzmann equation (0) can then be converted to diffusion 
equations in a standard way by doing a truncated moment expansion (see section 4). 

The largest contribution to baryogenesis in the MSSM comes from the chargino and neu- 
tralino sectors. For the charginos, the CP violating effects are due to the complex parameters 
1712 and in the mass term 

^.Mv. . Fo„ ( - to ( I I- 

The complex phases, combined with the mixing due to the Higgs fields, which vary inside 
the bubble wall, give rise to spatially varying effective phases for the mass eigenstates, which 
induce CP-violating currents for these excitations. To get analytic results, one can try to 
compute the current to leading order in an expansion in derivatives of the Higgs fields. This is 



the procedure followed in all methods designed to work on the thick wall limit [0, y, |T8 
This approximation cannot be used in the quantum refiection case 0, [1^, |^, |1^ , which can 
be relevant in the limit of very thin bubble walls. 
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We comment here on an apparent discrepancy in the hterature concerning the derivative 



expansion of the chargino source. References and obtained a source for the Hi — H2 
combination of Higgs currents of the form 

Sh,-h, ~ Im(m2 /i) {HiH!, - H,H[), (5) 



whereas ref. [18 found the other orthogonalhnear combination, ifiif2+-f^2-f^i- We previously 
beheved that the disagreement was because of fundamental differences between our CFM 
formalism and those of refs. [0, |17|, ^ . However we recently understood pO| , [2T[] that 



the difference was partially due to the fact that we were in fact computing the source for 
Hi + H2, for which the result is 

Shm ~ Im(m2 /x) {HiH'2 + H^H'i), (6) 

Therefore the disagreement about the sign was spurious: it can be shown that all three 
methods actually agree with eq. (P); it simply was not computed by the other references 



y,|i7|,^|i9i. 

The reason that the combination Hi + H2 was not considered by other authors is that 
it tends to be suppressed by Yukawa and helicity-flipping interactions from the term in 
the chargino mass matrix. Let us define chemical potentials for Hi, H2, left-handed third 
generation quarks ^3 and right-handed top quarks t, which we will assume are equal to 
the chemical potentials for the corresponding supersymmetric partners, as a consequence 
of supergauge interactions mediated by gauginos. If all the interactions arising from the 
Lagrangian 

V = iJ,hih2 + yh2URqL + yuRh2LqL + yu*Ji2LqL 

- yfihiqluR + yAtqLh2U*R + h.c., (7) 

were considered to be in thermal equilibrium, they would give rise to the constraints fiHi — 
/iQg + /iT = 0, + /^Qs — /iT = and fiHi + f^H2 = 0. If these conditions hold, the effect of 
the source Shi+H2 is clearly damped to zero. However, the rates of the processes coming from 
(0) are finite, and by studying the diffusion equations one can show that there are corrections 
of order (DhT)^^^'^ ~ 1, where we used Dh — 20/T and the Yukawa rate T ~ 0.02T (see eq. 
( |157|) and the discussion following). 
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Even in formalisms where the source Shi-H2 is nonvanishing [0, ^ |I9[, one should 
then not neglect the source Shi+H2 without first checking whether the other source, Shi-H2 
really gives a larger effect. In fact Shi-H2 does suffer from a severe suppression: quantitative 
studies of the electroweak bubbles in the MSSM show that the ratio H2/H1 remains nearly 



constant inside the bubble walls [^61 , |T2[; in Monte Carlo searches of the MSSM parameter 
space, the deviation from constancy is typically at the level of one part in 10^, and never 
more than 0.02. Therefore the source Shi-H2 is suppressed from the outset by a factor of 
10^ — 10^ relative to Shi+H2^ which is much worse than the Yukawa equilibrium suppression 
estimated above. 

In the CFM the situation concerning the source Shi-H2 is even worse: we will show that 
there will be no source arising from classical force, when computed correctly. To see this is 
actually quite subtle, and relates to the question of the gauge invariance we have referred 
to. If the problem is considered solely in terms of the canonical variables, there appears to 
be a non-trivial source of the form (H). That this term is unphysical however, is indicated 
by the fact that it can be transformed away by a field redefinition of the form 

Kl - e^"'/ii, (8) 

where ~ / lm{m2fi){H[H2 — H2Hi)dz. Below we will see that no such field redefinition has 
any effect on the physical momenta or currents, and hence should not give rise to a physical 
force (see also ||21|). In our treatment in terms of the gauge invariant kinetic variables this 



result is evident. In particular then the new source for baryogenesis in the CFM picture 



found in [^] is absent in our treatment. 

Our main result is that baryogenesis remains viable for a large part of the MSSM param- 
eter space, possibly with the explicit CP-violating phase as small as arg(m2/i) ~ 10^'^. The 
efficiency depends on the assumed squark spectrum, and the strongest baryoproduction cor- 
responds to the light right-handed stop scenario, which is also independently favored by the 
sphaleron wash-out constraint [l^, |Tl| . The resulting asymmetry has a complicated depen- 



dence on the wall velocity and for some parameters it peaks around the value of Vw — 0.01 
which has been indicated by recent studies of Vyr 



The rest of the paper is structured as follows. In section 2 we consider the simple case 
of a Dirac fermion with a complex spatially varying mass. We determine the dispersion 
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relation for the two helicities to leading order in Higgs field derivatives, and find from it the 
group velocity and the physical force acting on a fermion. We also compute and interpret 
the currents in the absence of collisions and show explicitly how the gauge invariant force 
can be identified from canonical equations of motion. In section 3 we employ the formalism 
in the case of the MSSM, in particular, computing the dispersion relations, group velocities 
and force terms for the charginos. (Squarks and neutralinos are also discussed here.) In 
section 4 we derive the diffusion equations, complete with the CP-odd source terms from the 
Boltzmann equations, using a truncated expansion in moments of the distribution functions. 
In section 5, these general results are applied to find and solve the appropriate set of diffusion 
equations which determine the chiral quark asymmetry in the MSSM. The rate of baryon 
production due to the excess of left-handed quarks is also computed in section 5, and our 
numerical results are given in section 6. In section 7 we present our conclusions, and a 
discussion of how the present results differ from previously published ones. 

2 Introductory example: Fermion with complex mass 

To understand some of the subtleties which arise when solving the equations of motion in 
the WKB approximation, let us first consider the example of a single Dirac fermion with a 
spatially varying, complex mass, 

(i7^9^ - tjiPr - m*PL)ip = 0; m = \m{z)\e'^^^\ (9) 

where Pr^l = |(1 =1=75) ^^e the chiral projection operators. We wish to solve eq. (P) approx- 
imately in an expansion in gradients of \m\ and 6. To simplify the solution we boost to the 
frame in which the momentum parallel to the wall is zero {p^ = Py = 0) and consider first 
positive energy eigenstates, ~ e"*"^*. Then, because spin is a good quantum number, we 
can write the spin eigenstate as a direct product of chirality and spin states 

= e-'"* [r])® = sXs, (10) 

where Rg and Lg are the relative amplitudes for right and left chirality, respectively (and 
we are using the chiral representation of the Dirac matrices). Spin s is related to helicity A 
by s = Xsign{pz). Inserting (^) into the Dirac equation (|^) then reduces to two coupled 
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complex equations for complex parameters Lg and Rg: 

{u — isdz)Ls = mRs (11) 
{u + isd,)Rg = m*Lg. (12) 

We can now use eq. (|l^) to eliminate Rg from (|T2D , which then becomes a single second 
order complex differential equation for Lg-. 

[uj + isdz) — {'^ — isdz) — m* \ Lg = 0. (13) 



m 

To facilitate the gradient expansion, we write the following WKB ansatz for Lg-. 

Lg=we'rp^('')d^'. (14) 

We have suppressed the spin index s in w and pc for simplicity. Inserting ([l^) into eq. (|13D 
we find the following two coupled equations (real and imaginary parts of (|l^)): 

9 , ,9 9 / N ^/ \m\' w' w" , , 

\m\ w w 

TfL ' 

2pcw' + p'w — j — -{suj + pc)w — 6'w' = (16) 
m 



While complicated in appearance, eqs. ([T^-|T6|) are easily solved iteratively. For example, to 
the lowest order one sets all derivative terms to zero, whereby the first equation immediately 
gives the usual dispersion relation cu^ = + |mp. It is also easy to extend the dispersion 
relation to first order in derivatives, because the contributions proportional to w' decouple 
from eq. (|15|) at this order: 

suJ+Po , , , . 

Pc = Po + Scp—^ + a , (17) 



where po = sign{p) u'^ — |mp. In(0) we have shown the generalization to antiparticles by 
including the sign Scp, which is +1 for the particle and —1 for the antiparticle. This follows 
from the fact that the Dirac equation for antiparticles is obtained from (P) by the substitution 
m —m*, which changes 6 to —9. The arbitrary function a'{z) reflects the ambiguity in 
the definition of momentum p^ in (|l^), because of the freedom to perform vector-like phase 
redefinitions of the field, — > e*"*-^-*-?/', which cause Pc Pc + ct'. This 'gauge' dependence 
reflects the fact that pc is not the physical momentum of the WKB-state, a quantity which 
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we will explicitly identify and show to be gauge independent below. In the preceeding, we 
considered the left-handed spinor Lg. The same procedure applied to Rs gives 

Pc=Po + Scp^ + a , (18) 



because of the sign difference between eqs. (|TT|) and (1121). The factors {su + Pc) are likewise 
replaced by su—pc in (plSj) and (|T6|) . In the following we will show that this difference actually 
does not have any physical effect: the group velocity and force acting on the particle is the 
same whether one uses (17) or ([181) . For simplicity we continue to refer to the relations for 



Ls unless the contrary is explicitly stated. 

2.1 Canonical equations of motion 

As anticipated above pc can be identified as the canonical momentum for the motion of the 
WKB wave-packets. To see this more clearly, let us first invert (|l7|) to obtain an expression 
for the invariant energy Q 

/ s9' 

^ = \/{Pc-acpy + \m\^-Scp — , (19) 

where a^p = a' + ScpO' /2 in the left- and acp = a' — s^pQ' j1 right chiral sector. (This 
difference in acp has no consequence what follows, which is why we have suppressed the 
indices referring to chirality). Identifying the velocity of the WKB particle with the group 
velocity of the wave-packet (corresponding to the stationary phase condition of the WKB- 
wave) it can be computed as 

/ON Vc~ 0(CP 

= (^Pc^)x = 



{Pc - OIcpY + |m|2 

Po( s\m\^e'\ 



^This discussion is closely analogous to the motion of a particle in an electromagnetic field, which can be 
described by a Hamiltonian 



R = v/(Pc - eA)2 + m2 + cAq. 



Here the canonical momentum Pc is related to the physical, kinetic momentum p = nvw j \J\ — = ujWg by 
the relation Pc = P + eA. Canonical momentum is clearly a gauge dependent, unphysical quantity, because 
the vector potential is gauge variant. Similarly canonical force acting on Pc is gauge dependent, but the 
gauge dependent parts cancel when one computes the physical force acting on kinetic momentum: 

Pfc = -d^H - edtA = e(E + V X B). 



8 



where the latter form follows on expanding to linear order in \m\'^9' /uj after eliminating 
Vc — C(cp with (^). Vg is clearly a physical quantity, independent of the ambiguity in 
definition of Pc- Given energy conservation along the trajectory we then have the equation 
of motion for the canonical momentum viz. 

Pc = -{dx^^)p. = Vga'^.p - - — ■ i07T + Sep— (21) 

[U-\-Scp — ) ^ 

which, like the canonical momentum itself, is manifestly a gauge dependent quantity, through 
the first term. Equations (^) and (|21| ) together are the canonical equations of motion 
defining the trajectories of our WKB particles in phase space. 

The physical kinetic momentum can now be defined as corresponding to the movement 
of a WKB-state along its world line 

P = OJVg. (22) 

This relation also defines the physical dispersion relation between the energy and kinetic 
momentum. We now calculate, using the canonical equations of motion (^) and (^Tf), the 
force acting on the particles defined as in eq. (|^) i.e. F = p = ujVg, where the latter follows 
trivially since ci; = along the particle trajectory. In particular we wish to verify explicitly 
that we obtain a gauge independent result for the force. Using the canonical equations of 
motion we have 

Vg = x{d^Vg)p^+ Pc{dp^Vg)x 

= Vg{d^Vg)p^ - {d^u)pXdp^Vg)^. (23) 

Using the form (pO]) for Vg, differentiating and substituting with the dispersion relation (|T 
we find 

.2 



{dxV. 



m 



gjpc 



(uj + Sapfr 



\m\\m\' 



from which it is easy to see that the gauge terms (in acp) cancel out exactly in (|23| ) and 
that the force is given by the gauge independent expression 

|m||m|'u; s9" \m\'^uj 
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CP 2 ) 



which to hnear order in 9' can be written as 



P = 



m\\m 



ijj 



CP 



2cj2 



(26) 



The force therefore contains two pieces. The first is a CP-conserving part, leading to hke 
deceleration of both particles and antiparticles because of the increase in the magnitude 
of the mass. The second part, proportional to the gradient of the complex phase of the 
mass term, is CP-violating, and causes opposite perturbations in particle and antiparticle 
densities. 

In connection with eq. ([TSD we mentioned the difference in definition of canonical momen- 
tum for left- and right-handed particles. From the immediately preceding discussion we can 
see that this difference gets absorbed into the definition of the unphysical phase acp- Indeed, 
for the right-handed fermions one should define a^p = a' — ScpO' /2 instead of a' + ScpO' /2. 
Since we have just shown that a^p cancels out of physical quantities, the difference between 
the dispersion relations derived from the spinors Lg and Rg has no physical effect. On the 
other hand, it is true that for relativistic particles will represent a particle with mostly 
negative helicity and Rg will correspond to a mostly positive helicity particle. The informa- 
tion about helicity (A) is contained in the spin factor, s = Asign(p^), and this does have a 
physical effect: particles with opposite spin feel opposite CP-violating forces. 

2.2 Currents 

Let us conclude this section by considering currents of WKB states under the influence of 
the CP- violating classical force (pG]). The current can be defined in the usual way. 



algebra one finds to this order the solution 

^p,, = I / \ I ® Xs e*/p«+4^^+-' (28) 

^2p+{uJ + Spo) \ \m\ V 'W^J J 

where = Po + slj9'/{2pq) and = + 9'/2. With this expression, it is simple to show 
by direct substitution that the current (|27|) corresponding to a WKB-state becomes 




(27) 



Now eq. (|T6D can be used to solve for w to first order in gradients. After some straightforward 




(29) 
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where we have restored the trivial dependence on p|| = {px,Py, 0) by boosting in the direction 
parallel to the bubble wall. This result confirms the intuitive WKB-particle interpretation; 
in the absence of collisions the quasiparticles follow their semiclassical paths, and when 
they slow down at some point the outcome is an increase of local density proportional to 
the inverse of the velocity. Because of this compensation of reduced velocity by increased 
density, the 3-D particle flux (j) remains unaffected by the classical force. 

Let us finally compute the current arising from a distribution of WKB-quasiparticle states 
using the physical dispersion relation. Under our basic assumption that the plasma can be 
well described by a collection of WKB-states we can write 

= J (2-) ^ (-^ - + ^-^) . (30) 

After integrating over pz, this becomes 

in perfect agreement with the result (^9|) . The current (|3ll) was recently derived from more 
fundamental principles in ref. |^ (see also |^); it was argued that the slightly more gen- 
eral result obtained in |^ would reduce to the form (|3lD in the limit of frequent decohering 
scatterings; this limit is of course an underlying assumption in the WKB quasiparticle ap- 
proximation used here. 



3 Application of WKB to the MSSM 

In this section we extend the previous analysis of dispersion relations and canonical equations 
to the case of the MSSM. The most natural candidate to effect baryogenesis in the MSSM 
would appear to be the left chiral quarks themselves, because any CP-odd perturbations in 
their distributions should directly bias the sphaleron interactions. However, in the MSSM 
the Higgs field potential is real at tree level, and therefore the CP-violating effect on quark 
masses arises only at one loop order. Moreover, the contribution from CP violation present 
in the supersymmetric version of the CKM matrix is potentially suppressed by the GIM 
mechanism, like in the case of the SM Excluding a direct source in quarks, one must 
look for CP-violating sources in various supersymmetric particles. These species include 
squarks, which couple to quarks via strong supergauge interactions, and charginos, which 
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couple strongly to the third family quarks via Yukawa interactions. We also comment on 
neutralinos, which have couplings similar to those of the charginos. 



3.1 Squarks 

After quarks the natural candidate to consider in the supersymmetric spectrum are the scalar 
partners of the third family quarks. The top squark mass matrix can be written as 

^^^^-[y{AH2 + fi*H,) ml ) ^^^^ 

in the basis of the left- and right-handed fields q = (tL, tR)^. Here the spatially varying 
VEV's Hi for the two Higgs fields are normalized such that in the zero temperature vacuum 
\/2{Hf + H2) = 246 GeV. The parameters uiq ^ refer to the sum of soft SUSY-breaking 
masses and VEV-dependent y^m^ and D-terms, but their explicit form will not be important 
here. Since squarks are bosons, they obey the Klein-Gordon equation 

{dl-dl + Ml)q = Q. (33) 

As in the case of a Dirac fermion, we first boost to the frame where the particle is moving 
orthogonal to the wall {p^ = Py = 0)- The chiral structure encountered in the fermionic case 
is missing here, but the problem is complicated by left-right flavor mixing. To deal with this 
mixing, at first order in the derivative expansion, it is easiest to perform a unitary rotation 
Uq to the eigenbasis of M|. The explicit form of the rotation matrix is 



Uq = diag(e^^-) , ,^ ^ ( ^(^^ + w/? A ^ h (34) 



A,(A, + A,)V -^l i(A, + A, 



where = y{A*H2 + fi*Hi), Ag = rriQ — mfj, and A = y A^ + 4|agP. The diagonal matrix 
diag(e*'^9*) contains arbitrary phases by which the local mass eigenstates can be multiplied, 
or equivalently the ambiguity in the choice of the rotation matrix, due to the U{1) gauge 
invariance of the lagrangian. After the rotation, eq. (^) becomes 

{00^ + dl - Ml + f/2 + 2f/i9,) qa = 0, (35) 

where Mj is a diagonal matrix, Ui = U^d^Ug and U2 = U^d'^Uq. We can formally write ( pSD 

as 

= 

= 0. (36) 
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The quantities -D±:p's are differential operators in tlie rotated basis, wliicli makes it impossible 
to exactly decouple the equations for the variables q±. However, one can show that they do 
decouple to first order in the gradient expansion. To this end we first write the equations 
(^) in the form 

(D++D__-Z}_+D+_)g_ + [D_+,D++]g+ = 

+ = 0. (37) 

It is easy to see that the commutator terms are of second order or higher in derivatives of mass 
matrix elements. Similarly, the products of the off-diagonal terms -D±^-D=p± are second order 
or higher and can be neglected. Finally, one can show that D^j^D±±q± = c± D±q± + 0{d'^), 
where c± are some constants. One then has simply 

D±±q± = (38) 

up to second order gradient corrections. Inserting the WKB ansatz q± = w±e^-f P"^'^^ into 
(p8|) one finds 

(2ipc±w'^ + ip'^^ + u;^ - m| - pj^ + 2ipc±Ui±±j w± = 0. (39) 
Breaking up the real and complex parts of the equations, we get 



uj''-ml-pl^ = 2pe±Im(f/i±±), (40) 



2pe±^ + ^Pe± = 2p,±Re([/i±±). (41) 
The correction term Ui±± appearing in the above equations is in fact purely imaginary: 

= ^ a^(a/+ ^y MA;f,){H[H, - H^H,) + , (42) 

where (p'^^ = (f/gdiag(0^j)f/j)-|-± are still some arbitrary phases. Using this notation, we 
see that the dispersion relation acquires an energy-independent shift from the leading order 
result: 

Pc±=Po±-^;±, (43) 



where po± = sign{pz)y uj"^ — m\. Curiously, the parametric form of the shift (|42|) is the same 
as what appears in the source derived for squarks in |2^. In our method this correction 

13 



originated from a local rotation in the flavour basis of the mass eigenstates. Similarly, in 



25| the source oc lm{Alfi){H[H2 — H2H1) was found by performing an expansion to a finite 
order in the flavour nondiagonal mass insertion over temperature, which is an approximate 
way of taking into account a rotation of eigenstates in a varying background. In the present 
context we can see, however, that this shift is unphysical, because of the arbitrariness of 
the phases (p'^^. in (0). For example, they could be chosen to make 6'^_|_ = 0. This is only 
possible because the expression ( ^21) is a function of x only, and not p. Indeed, proceeding 
in analogy with the fermionic case of section (2.1), we find that pc± is to be identified as the 
canonical momentum of the system. Moreover, defining the physical momentum p± through 
the group velocity as in (p2D, we find 

p± = ujVg± = po±. (44) 

Thus the physical momentum gets no corrections to first order accuracy. This implies that 
neither does any classical force arise at first order in gradients. Notice also that, because 
Re(f/i-i-±) = 0, the normalization of the state can be computed only to the zeroth order from 
(PI), which gives w± = C±/y/po±, where C± are constants. 

Because the CP- violating source can only arise at second order in the gradient expansion 
in the squark sector, it is parametrically small compared to a fermionic source (to be derived 
for charginos below). Given the range of wall widths compatible with a sufficiently strongly 
first order phase transition in the MSSM [^, |l^, we can estimate this suppression roughly 
to be of the order ~ km' /m ~ 1/3T£^ ~ 1/30 for a particle with thermal de Broglie wave 
number ~ 1/3T and wall width 1^ ~ 10/T. We will accordingly neglect the squark source 
henceforth. 

3.2 Charginos 

An asymmetry in Higgsinos is efficiently transported to left-handed quarks via strong Yukawa 
interactions with third family quarks. The chargino mass term, 

^ijM^i + h.c, (45) 

contains complex phases required for a CP-violating force term. In the basis of Winos and 
Higgsinos the chiral fields are 
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= {Wt.hlLf (46) 

and the mass matrix is 

where the spatially varying VEV's Hi are definded as in the squark case. The corresponding 
Dirac equation, in the frame where = = 0, is 

{il^dt - i-fsd, - M^Pl - MPn)^ = 0. (48) 

To solve it in the WKB approximation, we follow the same procedure as with the single 
Dirac fermion and the squark cases above. First we introduce the spin eigenstate as a direct 
product of chirality and spin states, where spin s and helicity A are related by s = A sign{pz): 

= e-*"* [r])^ ^3$. = s $ (49) 

In contrast with the simple Dirac fermion, the relative amplitudes of left and right chirality, 
Ls and Rs, are now two-dimensional complex vectors in the Wino-Higgsino flavor space. 
Keeping this generalization in mind, the solution proceeds formally in analogy to the case 



of a single Dirac fermion; inserting (^91) into the Dirac equation (48) gives 

{uj-isd,)Ls = MR, (50) 
{uj + isd,)Rs = M^Ls. (51) 

From eq. ( ^0] ) we have Rg = M^^{uj — isdz)Ls, which when substituted into (|5T| ) gives 

(cu^ + - MM^ + is{Md,M-^){uj - isd,)) L, = 0. (52) 

Since Lg is a two-component object, writing the WKB-ansatz is somewhat more involved 
than it is for a single fermion. But since MM'^ is a hermitian matrix, we can rotate to its 
diagonal basis, similarly to the squarks. Eq. ( |52D then becomes 

(io^ + dl -ml + 2Uid, + U2 + tsAiioo - isd,) + A2) = 0, (53) 

where the superscript in indicates that we are in the basis where MM^ is locally diagonal 
as a function of distance z from the wall. The 2x2 matrices (in the Wino-Higgsino flavor 
space) in (|53D are defined by 

Ai = U{Md,M-^)U^; A2 = AiUi. (54) 
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The explicit form of the rotation matrix U which diagonahzes MM^ is similar to the one 
encountered in the squark case: 



^ /A(A + A) V i(A + A) y' ^ ' 



where 



a = m2Ui + /i*'U2 

A I |2 I |2 I 2 2 

A = \m2\ — + ^2 — 1*1 



A = ^A2 + 4|a|2 

Ui = gHi. (56) 

The arbitrary angles 0j will eventually enter the dispersion relation as physically irrelevant 
shifts in the canonical momenta, similarly to the squark case and the case of the single Dirac 
fermion. The diagonalized MM^ matrix, mj^ = diag(m^,m^), has the eigenvalues 

"^l = ^(l"^2|'+|/u|' + «2 + w?)±^- (57) 

Broken into components, labeled by ±, and suppressing the spin index on , equation 
(p3|) can be written as 

D^^Lt- D^+L% = 

D++L%- D+_Lt = 0. (58) 

Just as in the squark case, one can show that the mixing terms in (|58|) can be neglected to 
the first order in the gradient expansion, and it is sufficient to solve the decoupled equations 

D^^L'^ = 0. (59) 

Inserting the WKB ansatz, Lj_ = w±e*^^='='^^, into (^), and writing -D=p^ explicitly, we obtain 

(oo'^ - P±- m\ + zp'± + '^W±dz + 2zp±f/i±± + is{uj + sp±)AiJ^j^ w± = 0. (60) 

Taking the real and imaginary parts of this equation we have 



cu^ - p| - m| = Im (2p±f/i_t± + s{uj + sp±)Ai^j.) (61) 



p'± + ^P±— = Re (2p±f/i±± + s(c^ + sp±)Ai±±) . (62) 
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Equations (|6ll - |62|) are similar to the equations ( ^0[]4l|) for squarks apart from the appearance 



of new contributions from the matrix Ai in (|6l| - |62D ; as we shall see, this difference is crucial. 
Eq. (|6l|) gives the dispersion relation to first order in gradients; however (^) gives the 
normalization (wj-) only to zeroth order, because integrating w'_^ eliminates one derivative. 
To this order w± give the usual spinor normalization, but with a spatially varying mass terms. 
If we needed to know w± also at first order, it would be necessary to include second order 
corrections to (|62D. Luckily we do not need these results here, and therefore concentrate on 
the dispersion relation in the following. 

To find the leading correction to the dispersion relation we need to compute the diagonal 
elements of the matrices Ui and Ai. 

T A Im(m2/i) , , 

21m ( 1712 fi) , t i\ ~t 

Im?7i±± = ± ^ {U1U2 - M2M1) + ^0L±^ (63) 



where ^'^^^ = (t/diag(09?7''')±±. Defining po = sign(p)-y'u;2 — |m±P the dispersion relation 
for the states associated with becomes 

= P0± T ^CP^^^ -^{U,U2 + U2U,) 

T Sap^^^^{u,u'2-U2u[) ± ^A4>',^. (64) 

The sign Sijp is 1 ( — 1) for particles (antiparticles). The signs i refer to the mass eigenstates; 
below, we will want to focus on the state which smoothly evolves into a pure Higgsino in the 
unbroken phase in front of the wall. This will depend on the hierarchy of the diagonal terms 
in the chargino mass matrix in the following way: 





1^1 


> 




1 


1^1 


< 


■^2 



Ph. = y r : L (65) 



The reason for identifying p with the Higgsino /i2 is that it plays the role of the left-handed 
species in the mass term, as we have written it in eqs. (^5| - ^6D . 

In the diffusion equations to be derived in the following sections, we will treat the 
charginos as relativistic particles, whose chirality is approximately conserved. We there- 
fore would also like to know the dispersion relation for the other flavor component, h2, which 
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is associated with the right-handed spinor i?f . Our convention for the sign of /i2-number is 
the supersymmetric one, where the Higgsino mass term is written in terms of left-handed 
fields only and has the form 

fihih2 + h.c. (66) 

Explicitly, we identify 

hi - htL = ChUr 

h2 ^ /l2,L (67) 

That is, hi is identified with the CP conjugate of hf j^, whereas hf j^ itself is the particle 
represented by the spinor R^. Therefore we must remember to perform a CP conjugation of 
the i?f -field dispersion relation if it is to represent the states which we call hi, a point which 
can be somewhat confusing.^ 

Going through the same steps as for to find the dispersion relation for i?f , we obtain 

u;^ - - m| = Im {2p±Vi±± + s{uj - spk)Bi^±) (68) 

where Vi and Bi respectively can be obtained from Ui and Ai by exchanging ui ^ U2 and 
taking complex conjugates of m2 and /i. Taking into account the additional sign change 
Sep — >■ —Sep required by our convention for the meaning of hi, we obtain 



R siuo — spQ±)\m{'m2U) , , 

P± = P0± T Sep — -^—{uiU^ + U2U,, 



± Sap^^^{ulu'2-U2u',)T^A4>'Hi. (69) 

where A = |m2p — + u\ — u\, = (Vdiag(09V^)±±, and by definition Scp = 1 for 
the state hi^i. Similarly to (|65|), we identify 





1^1 


> 


■^2 


\ pt 




< 





P'hi { „R |,,| ^ |rr)„l ^^^^ 



The term proportional to Uiu'2 — U2u'i in (|6^) and (|69|) was not included in our earlier 
work. (The complete dispersion relation was however given in reference recently.) This 



term is odd under exchange of the higgs fields, and appears potentially viable to produce a 
source Shi-H2 i^i the diffusion equations for the combination hi — h2. Indeed, if one does not 
•^We erred on this point in 



keep in mind that the momenta p^'^ in (^) and (|69D are the canonical momenta, and not 
the physical momenta, one is easily led to infer (as recently in ref. [^) that there is a CP- 
violating force, since the canonical equation of motion = —{dxUj)p^ includes a contribution 
from this UiU2 — U2u[ piece. As discussed in section 2, however, the latter quantity is, like the 
canonical momentum itself, a gauge invariant quantity which changes under arbitrary overall 
local phase transformations on the fields. Just as in the squark case, the U1U2 — n2^^'i-part of 
the dispersion relation is energy- independent, can be absorbed into the arbitrary phase factor 
A(p'j_, and as such does not represent a physical quantity. In fact, grouping all unphysical 
constant terms from the r.h.s. of ( |64D or ( |69D into a common arbitrary phase factor, we get 

Pi'"* = Po± + ScP^^ + a'^", (71) 
with the physical phase 6'_^ defined by 

e'^^T'^^^^M + u,u[). (72) 

This result shows how the chiral force depends on having both CP-violating couplings in 
the Lagrangian (in this case the phase of 1712 fi) and spatially varying VEVs — otherwise 
the phases could be removed by global field redefinitions — so that the force is operative 
only within the wall. Treating the whole problem in the kinetic variables, as we do here, 
one avoids by construction the problems of gauge variance one encounters when using the 
canonical variables. 

The physical part of the dispersion relation is identical for both species of higgsinos, hi 
and h2- Hence these states will have the same group velocities and experience the same 
physical force in the region of the wall. The form ( [7TD for the dispersion relation is also 
identical to that for the single Dirac fermion, with the simple replacement ^ — 6'-|-, so that 
we immediately have 

Po± sscpm\e'^ 

- -^ + ^K«i)'. (74) 

Since the force terms are identical for both kinds of higgsinos, so will be the source terms 
in their respective diffusion equations. This will become explicit when we prove in the next 
section that the source is proportional to a weighted thermal average of the force term. 
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The outcome is that the hnear combination Shi-H2 considered in and in |Tj, is not 
sourced at all in the classical force mechanism, at leading order in the WKB expansion. 



3.3 Neutralinos 

Neutralinos are an obvious candidate to study after charginos. The mass term for neutralinos 
can be written analogously to that of the charginos as ^Mn"^ l+ h.c, where in the basis of 
gauginos and neutral Higgsinos, = {B, W^, j^, hj^j^'^, 

V = mz\ ■ a ■ a -a /? ('5) 
\ — sm tly^ sm p sm O^j cos p / 

where tan/? = V2/V1, and mz is the Higgs- field-dependent Z boson mass. Because it is a 
4x4 matrix it is more difficult to solve for the WKB eigenstates of the neutralinos than for 
the charginos. However the structure of the mass matrices is sufficiently similar to suggest 
that the chiral force on neutralinos and hence the magnitude of the produced asymmetry in 
neutralinos is not in any way parametrically different from that of the charginos, and should 
be quantitatively similar as well. However, the transport of the asymmetry from neutralinos 
to the quark sector is much less efficient due to smaller gauge-strength coupling to fermions. 
We will therefore limit ourselves to a computation of the chargino contributions alone in the 
following estimate of the baryon asymmetry. 



4 Transport Equations 

Our next task is to determine how the nontrivial dispersion relations lead to CP-odd pertur- 
bations on and around the bubble wall. In particular we need to determine how an asym- 
metry in left-handed quarks is produced which drives the electroweak sphalerons to generate 
the baryon asymmetry. Indeed, our primary source particles with direct CP-violating inter- 
actions with the wall (charginos) experience no baryon number violating interactions, and 
therefore the CP-violating effects must be communicated to the left-handed quark sector via 
interactions. Within the WKB approximation the plasma is described by a set of Boltzmann 
equations for the quasiparticle distribution functons. We will not attempt to solve the full 
momentum dependent equations; instead we will use them as a starting point to derive, by 



20 



means of a truncated moment expansion, a set of diffusion equations for the local chemical 
potentials of the relevant particle species. 

The advancing phase transition front (bubble wall) distorts the plasma away from the 
equilibrium distributions which would exist if the wall were stationary. The exact form of 
the distortion is complicated, but a simple ansatz can be made in the present situation for 
two reasons. First, the perturbations in the chemical compositions are small in amplitude, 
because they are suppressed by the presumably small phase in the CP violating part of the 
force exerted by the wall on the particles. Second, the elastic interactions enforcing the 
kinetic equlibrium are much faster than the ones bringing about the chemical equilibrium 
and moreover, they are also very fast compared to the wall passage time-scale. 

We first need to determine the form of the local equilibrium distribution function for 
the WKB-states. To this end we need to more accurately specify what we mean by the 
particle interpretation of WKB-solutions, i.e. what is the appropriate local ^-component of 
the momentum. Following our reasoning in the treatment of the flow term, we argue that for 
an interaction with a mean free path less than the wall width, F^^ the WKB-state 

can be approximated by a plane wave (particle) solution 

V'WKB^ exp[i{ujt - p|| ■ X|| - p^z)], (76) 

where p|| is the momentum parallel to the wall and pz is the kinetic momentum defined 
by the physical dispersion relation such as (pO|) or (ffSl). If F~^ << one can to a good 
approximation compute the collision terms by extending this wave solution to infinity, which 
then formally implies the usual Feynman rules for the states with conserved 4-momentum 
[uj; p\\,Pz)- Under these considerations the appropriate equilibrium solution in the wall frame 
is given by 

^-°(P'^) = e/37.(^+^^p.)±r ^^^^ 
where v^, is the velocity of the bubble wall, 7^, = l/^Jl — t>^, (3 = 1/T and p is the kinetic 
momentum. The signs ± refer to fermions or bosons, respectively. It is clear that the 
distribution ( [77D makes the collision integral vanish for the states ([76| ) with a constant pz. 
However, if F~^ is not small in comparison with the wall width, one should expect corrections 
of the order ~ Tdz to the collision terms. These corrections are parametrically of the form 
of the spontaneous baryogenesis terms discussed earlier in the literature |15], ^ |l8l. We 
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do not attempt to derive these terms here, because based on our earher estimates [|T8l, we 
expect them to be subleading to the source term deriving from the classical force in the flow 
term. 

The need to define the equilibrium distribution externally is an inherent shortcoming 
of the semiclassical method. In previous papers j^, a slightly different approach was 
taken, where the ansatz was parametrized in terms of the canonical momentum. In canonical 
parametrization however, one is faced with a more subtle task to obtain phase reparametriza- 
tion invariant physical results. We will give a detailed comparision of the two approaches in 
appendix B; here we comment that our present definition has a more direct physical interpre- 
tation and that the difference between the two approaches would lead to only small changes 
in our final results. At any rate, a more complete derivation of the transport equations. 



using the methods developed in [^, will be needed to settle the issue self-consistently. 
We now extend the ansatz ( [77|) by allowing perturbations around the equilibrium: 

where the local chemical potential function fii{x) and the local momentum-dependent devia- 
tion function 6fi{p, x) are expected to be small in magnitude. The first part of ( [78| ) describes 
the system in kinetic, but possibly out of chemical equilibrium, whereas the momentum- 
dependent part 6fi{p,x) describes any deviations from kinetic equilibrium. By definition of 
the local chemical potential f^i{x), the total number density is given by the first term of ( [78|) 
alone, and hence 

Jd^p6f,{p,x) = 0. (79) 

The 5f{x,p) term will play the role of an auxiliary field in the following derivation. It cannot 

be neglected, for the Boltzmann equations with a force term can only be self-consistently 

solved if there is a perturbation which is anisotropic in momentum space.0 The shape 

of 6fi{p,x) cannot be consistently restricted beyond (|79|) , since this would require some 

hierarchy between the elastic interaction rates that would allow certain perturbations to be 

damped more slowly than others. 

We now return to our derivation of the transport equations, adopting the ansatz ([78|) . 

We have derived the dispersion relation in the rest frame of the bubble wall, where it is also 

^The same is true of the derivation of the diffusion current j — —DWn in the simplest free diffusion 
approximation, or Ohm's Law in the presence of an electric potential. 
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easiest to derive the flow term of tlie Boltzmann equation. Tlie collision term is simple in 
the rest frame of the plasma however, and to exploit this we must Lorentz-transform the 
flow term to the plasma frame in the end. Since the wall velocities are expected to be less 



than 0.1 |28|, we will expand in and keep terms to first order in v^; hence 7^„ = 1 and 
it is sufficient to perform a Galilean transformation for the flow term. 

In the following subsection we will derive the diffusion equations in detail. The result 
will be a set of coupled equations for = /ij/T of the form — -Dj^f — v^^'^ + = Si, 
where Di is some yet to be determined diffusion coefficient, Tfj is a matrix of rates at which 
particles are lost due to decay and inelastic collision processes, and Si is the source term 
due to the chiral force. Application of this formalism to the MSSM will be considered in the 
next section. 

4.1 Derivation of diffusion equations 

Moment expansion. Our starting point is a set of coupled Boltzmann equations for the 
relevant particle species. Since the expansion of the universe is negligible compared to the 
diffusion and interaction time scales, we have simply 

(dt + y^g-d^ + F- dp)fi = C[f„ /„ ...]., (80) 

where is the group velocity and F = p is the classical force, as derived in sections 2 and 3, 
and the partial derivative with respect to x is now taken at fixed kinetic momentum p. Q To 
derive the diffusion equation for /ij, we first insert the ansatz (^) into (|80|) , approximating 
7^ = 1, which gives 
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[(d,uj),^-fi')^ + dJf,]+F,,,[{idp^uj), + v^)^ + dpM] =C[/„ /„...]. (82) 



where primes mean d/dz and the expression dfi/du refers only to the non-(5/j part of /j. It 
should be noted that the partial derivatives with respect to Pz in (^) are taken with respect 
to the kinetic momentum, and therefore do not correspond to the usual canonical equations 
^The Boltzmann equation was earlier ^] written in terms of canonical variables: 

(dt + vg • a, + • dpjf, - c[/„ /„ ...]., (81) 



Equations (pO|) and (81) are actually equivalent, since Eqn. (pOl) can be obtained from (pl\j by a simple 



change of variables {x,pc) {x,p{pc,x)). 
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of motion, i.e. {dp^uj)^ 7^ Vg and ~{dzUj)p^ 7^ Fj^^. However, these terms still cancel because 
their sum is a total derivative of u: 

Vg{dzu)p^ + Fi^z{dp^uj)z = z{dzuj)p^+pz{dp^uj)z = du/dt = 0. (83) 

Moreover, because the correction 5f is generated by the force, one can easily see that the 
CP-even correction arising from the Fi^zdp^Sfi term is of second order in gradients, and the 
CP-odd correction of third order, both one order higher than corresponding terms coming 
from the leading v^Fi^^ term. We can thus neglect it and we are left with 

^ [v^F,,. - /i'l) + ^5f: = C[f., /„ ...], (84) 

where we also wrote Vg = Pz/^^, which is an exact relation in the physical variables. As 
expected, the only source termQ in (p^ ) is proportional to the wall velocity, so that the 
trivial unperturbed Fermi-Dirac or Bose-Einstein distribution, /o = (e^*^ ± 1)^^ becomes 
a solution in the limit of a wall at rest in the plasma frame, as it should. We now wish 
to transform to the plasma frame, since the collision integral takes a simple form there. 
Since the wall velocity is small []27| , |28[| , we can make a Galilean transformation of (|8^) , 
Vg ^ Vg + Vw We can also replace fi by the unperturbed distribution /q to leading order in 
the perturbation. Then 

?7 + ("^^ ^ + ^ ^ ^^^^ 



for the Boltzmann equation (|80|) in the plasma rest frame. The term in (^) containing 
the force Fi is what pushes the distributions out of both kinetic and chemical equilibrium. 
Notice that the chemical potential in this equation is not a purely CP-odd quantity. It is 
actually a j»seu(iochemical potential [^] which also has a component caused by the CP-even 



part of the classical force. The latter is relevant for the dynamics of the wall expansion, 
and determines the wall velocity and shape |3y, |2^, For baryogenesis however, only the 



CP- violating part is important. We will therefore neglect the CP-even parts of fi and Fi in 
what follows, denoting by 6Fiz the CP-odd part of the force. 



^i.e. the term which is nonvanishing when ^ = Sf = 0. Notice that the coUision term vanishes in this 
hmit. 
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We next integrate both sides of ( ^5]) over momentum, weighting by 1 and by Pz/^-, 
respectively, to obtain two independent moment equations:^] 

where the average over momenta is defined as 

- 1 f Sd'pX'^. X=vll^\ 
^ '~T\dH^^\ld^vX. X = allothers. ^ ' 

In deriving (^6]) we used the fact that some terms are total derivatives of momentum, which 
integrate to zero. To leading order in the WKB approximation, we can take {{pz/u!)6f-) ~ 
{{Pz/^)Sfi)' since the error is third order in derivatives for CP-odd quantities, ~ m'9",m"9' 
(fourth order in the diffusion equation). This lets us express the l.h.s. of ( |5BD in terms of two 
variables 



= fii/T; 

v,,z ^ [^Sf^y (88) 

Writing further Pz/uj = Vp^, (the z-component of the physical group velocity) eq. (^) takes 
a more transparent form which shows its dependence on powers of the various velocities: 

-Vn^C + Kz = (Ci) (89) 
-{vlX + vX^, + v^{vpJFjT) = {vp^a) (90) 



To proceed, we express the moments of the collision integral appearing in equations ([89- 
Op in terms of the variables C,i and Vi^z- This is possible because all perturbations have small 
amplitude and we can linearize the collision integrals in fii and 6fi{p,x). In the Appendix it 
is shown that 

{vp^Q) ~ v,,zTl; (91) 

(Q) ^ -rf^Eef. (92) 



^This choice is not unique. To get the most reUable truncated equations, the moments should be chosen 
such that they sensitively probe the momentum region that is most relevant for the problem. Here the 
physical effect is not confined to any sharply cut region in the momentum space, and since we are performing 
an expansion in the wall velocity, it is natural to use the lowest moments of the unperturbed particle velocity. 
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where Tj is the thermally averaged total interaction rate for the state i, Tfj^ is the averaged 
interaction rate corresponding to an inelastic reaction channel labeled by k (we assume a 
sum over all allowed interactions) and the abbreviated notation J2j ^j'^ in (EU) represents the 
signed sum over the chemical potentials of particles participating in that reaction.^ 



The first moment, eq. (91) is proportional to the velocity (kinetic) perturbation. As 
expected, the relaxation scale for these perturbations is set by the total interaction rate. 
The zeroth moment, on the other hand, only gets contributions from inelastic interactions, 
because (Cd) = 0. The damping term appearing in (Cj) causes relaxation towards chemical 
equilibrium, which is reached when the signed sum of the chemical potentials, J2j goes 
to zero. The relaxation scale is set by the inelastic interaction rates F^, and is therefore 
much longer than that of kinetic perturbations. These terms also provide the couplings of 
particle species which leads to evolution of the asymmetries from the source species into the 
left-handed quark fields, a crucial element of the present baryogenesis mechanism. 



Diffusion equation. Differentiating the first moment equation (|9(]|) once, using the zeroth 
moment equation (|89|) to eliminate Vi^z in the resulting equation, and neglecting terms of 
order C(f^) (since v"^ <ti {VpJ) one obtains the following second order equation for ^j: 

- - v^C - E Y^^f' + rf.E = -"-^i^pMY- (93) 

Because of the normalization chosen in (|87|) , (see also (|171D ), an interaction coupling certain 
species may appear with different effective strengths in different coupled equations, depend- 
ing on whether the equation describes a fermion or a boson chemical potential, and what 
the particle's mass is. For later purposes it will be more convenient for a given interaction 
to have the same rate in all equations. We therefore multiply each equation by a factor 



where dp = and ± is + if particle i is a fermion and — if a boson. Each interaction rate 
is then rescaled using 

f ^ K,Ti. (95) 

^For example a term induced by the decay process i-^jk is more precisely written as Ti—,jk{^i ~ + Cfe)- 



26 



The rescaled decay rate is equal for a given reaction, regardless of which particle's equation it 
is appearing in, and the information about the differences in statistics and masses is encoded 
in Kj's. (In the massless limit, = 1 for fermions and Ki = 2 for bosons.) Also, as explained 
above, T^^/F* is small so that we can drop the third term in (|93D . We then obtain the usual 
form of the diffusion equation 

- «:,(AC" + Vv.O + E = (96) 

j 

where the diffusion coefficient Di, related to the inverse of the total scattering rate, is given 
by the usual expression: 

r* 3 ^ ' 

S.^-K.^{v,M)'- (98) 



and the source is 



The proportionality of the force term to the diffusion constant reflects the fact that the 
particles must be able to move in response to the force in order for it to have an effect. The 
larger the diffusion constant, the bigger is the response. 

5 Transport in the MSSM 

We now apply the general formalism given above to the MSSM. The goal is to see how 

the asymmetry in source particles due to the chiral force gets transmitted into left-handed 

quarks, which drive the sphaleron interactions to produce a baryon asymmetry. The network 

of diffusion equations connecting the various species of particles in the thermal bath can be 

simplified considerably if we take account of the hierarchy of reaction rates for inelastic 

processes that change particle identities. We take the electroweak sphaleron rate, of order 

a^T, to be the slowest of all the relevant interactions, hence the one we treat last (section 6), 

whereas the gauge interaction rates (order a^T) are fast and can be taken to be in equilibrium 

on the time scale -D/f^ for particles to diffuse in front of the bubble wall.0 

Intermediate between these extremes are certain inelastic processes that could possibly 

drive the chiral asymmetry and hence the baryon asymmetry to zero if they were sufficiently 

^ The estimate D/v^ can be obtained by setting the distance a particle diffuses in time t equal to the 
distance the wall translates it the same time i.e. y/Dt ~ v,,,t. We will sec that this is indeed the charatccristic 



scale in the Green's function (137) we derive below to describe the diffusion. 
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fast. These include strong sphalerons (F ss ~ ^t'^) s^iid helicity-flipping interactions T/j/ ~ 
(m^/T^)rt, the first of which tends to erase the chiral asymmetry of the quarks while the 
second tends to damp the chiral asymmetry of the charginos which is originally produced by 
the classical force. 

The only source for the asymmetry in our equations comes from the chargino sector. It is 
clear from our results in section 3, that the wino-like and higgsino-like mass eigenstates have 
equal and opposite sources. However, whereas the asymmetry coming from the higgsino- 
like states can be efficiently transported to a chiral asymmetry in quarks and squarks via 
strong Yukawa interactions, the asymmetry from winos can only transport via mixing effects 
through these same interactions. Indeed, the dominant gauge interactions between winos 
and quarks do not give rise to a transport of asymmetry, but produce a large damping 
term for the wino asymmetry. Moreover, because the left- and right-handed Winos have 
opposite source terms, it is only the difference of their chemical potentials, —iJ,^^+, whose 
diffusion equation has a nonvanishing source; and the thermal equilibrium of the supergauge 
interactions enforces the constraints //^+ — oc //^ oc //^s oc /i^s- Because the neutral 
gauginos have helicity-flipping Majorana masses which tend to drive their chemical potentials 
to zero, this brings another, gauge-strength suppression to the wino chemical potential. 

For the sake of tractabihty of the problem, we will here neglect the chargino mixing 
effects in the interaction terms, which allows us to drop the winos from our diffusion equation 
network. Based on the above discussion, we believe this to be a very good approximation 
everywhere except in the region of parameters where mixing in expected to be large i.e. , 
when \ni2\ — |/i| ^ gHi. Even in this case, however, we err only within the bubble wall 
where the mixing is large; in front of the wall where (Hi) = 0, there is no mixing. Since the 
diffusion length of the Higgsinos is significantly greater than the wall width, (see Fig. 1(b)) 
the small mixing approximation is a good one for most of the region where has support. 

Another simplification we have made is to assume that the gaugino helicity-flip inter- 
actions are in equilibrium; this leads to the constraint that all particle chemical potentials 
are equal to the corresponding ones for their superpartners. Although neither of these two 
simplifying assumptions are particularly good approximations, we expect them to introduce 
multiplicative errors of at most of order unity in our estimate of the baryon asymmetry, 
without changing the important parametric dependences. 
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The relevant diffusion equations are most easily written in terms of the rescaled (s) quark 
and Higgs(ino) chemical potentials, C,i = fJ'i/T. Because of our assumption of vanishing 
gluino chemical potential and equilibrium of supergauge interactions, C,q = C,q for each light 
squark species. Moreover, the chemical potentials of quarks within a doublet are taken to 
be equal due to weak gauge interactions, for example ^f,^ = = ^gs- The asymmetries 
between the particle and antiparticle densities, plus those of their superpartners, include the 
Kj factors defined in (0): 

- ^Q3 = Y (2 + t^bL + '^h) ^93 . (99) 

and likewise for the other species. The factor 2 in ( p9t ) comes from the quarks, while the 
K factors count the squarks. (Recall that Kg = 2 if the squark q is light enough to be in 
equilibrium in the plasma, and Kg = otherwise.) Defining the diffusion operator 

A = -6(^A^+t^«,^j (100) 
the network of diffusion equations can be written as 

i^h^h + r.^fe.-e.3+6J + r./(4,+4J = Sh, (loi) 

^h^h, + (r.+r,^)(a,+e,3-6J + r./(4,+4J = Sh, (102) 
l{2+K~,^+KijVg^g, + (r,+r,^)(4^+e,3-6J -r,^(e/.,-e,3+6J 

- T^e{-x) (6« - e.a) + 2r,, = (103) 

+ Tme{-x) (6, - - r,, = (104) 

l{2 + K^^ + ^^jVgCg,,, + 2r,, = (105) 

1{i+kMU- r.. = o, (106) 

where stands for any of the right-handed (s)quarks except for the top. In section 3.2 we 
showed that both species of Higgsinos experience equal chiral forces, so that Shi = Sh2', we 
will denote it by Sh henceforth. The rates Fj are associated with terms in the interaction 
Lagrangian (|^: 

Thf ^ /i/ll/l2 
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Ty ^ h2URqL + yuRh2LqL + yuyti2LqL 

TyA ^ yAtqihun (107) 

In addition to these, Tss is the rate of strong sphaleron interactions, which interconvert left- 
and right-handed quarks of each flavor; we have used the shorthand 

Tss = Tss{2^q^ + 2^q2 + 2^q^ - " ^(Ir " ^cr " ^sr " ^bR " 6^)5 (108) 

to simphfy the appearance of the equations. Also TmO{—z) is the spin-flip rate due to the 
quark mass in the broken phase behind the bubble wall, which we take to be the region z < 0. 
{6{—z) is the step function.) There exist similar gauge-strength interaction terms damping 
the chemical potentials in the Higgs condensate region. They could easily be included in 
the analysis, but because they are subleading in strength, we omit them for clarity. 

Let us now introduce a "reference chemical potential," ^q, which satisfies the equation 

V,^Q - 2r,, = 0. (109) 

Then from ( [L05| ) and ( |106| ), one can express ^^^ j and in terms of ^q: 

U = Tx;— (110) 

i -h Kq^ 



Furthermore, adding ( |103D and ( |104D we get the equation 



Dq {^tR + 1(2 + + ^ijQ + 2Tss = 0. (Ill) 
This is again of the form (|109|) and can be used to infer the constraint 

^tR+2Rs^q, = ~^Q. (112) 

where we have defined the frequently appearing quantity 

R,^l(2 + t,-,^ + ^-J. (113) 



The condition ( |112|) is in fact equivalent to requiring the vanishing of total baryon number 



-Btot oc I]g(l + tiq)^q = 0, which is built into our equations until we explicitly introduce 
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the B-violating electroweak sphaleron interactions. These constraints leave us with four 
undetermined potentials, say ^t^, ^ga, ^h2- '^^^ further reduce this number by 

assuming the strong sphaleron interactions are in thermal equilibrium. Using ( |110| ) and 
( |112|) , the strong sphaleron term (|108|) can be expressed in terms of just and 



r = r 

^ ss ss 



2(1 + ksa)^g, + (a - 1)6, (114) 



where 

a=Yl + E (115) 

i=l,2 2 + f^u^L + l^diL q^t 1 + '^^R 

Assuming strong sphaleron equilibrium allows us to set ( |114| ) to zero, and thereby determine 

6. = 2^^^^e.3 = (l-C3)e.3- (116) 
Eq. ( |116| ) must be used with care when eliminating 6, from our system of equations. If 



used indiscriminately in ( |103| ) or ( |104| ) individually, two incompatible equations would result. 
The correct procedure is to form the linear combination of ( p.03| ) and ( |104| ) which does not 
contain the large term r^^. Eliminating 6, from this combination gives the correct result 
at leading order in an expansion in l/r<js. In this way we reduce our system of diffusion 
equations to a set of just three: 

T^hiK, + r,,(4,-c3e,3) + r./(4,+4J = ^H (117) 

T^hi-K, + (r. + r,^)(4. + c3e.3) + r,/(4,+aj = 5// (ns) 

a%3 - (rj/ + ryA)(42 + ^3^©) 

+ r,^(e/,, - c,i,,) - C3 r^e,3 0{-z) = 0, (119) 

where a = 1 — + k^. Since only the linear combination + is sourced, it is useful to 
define the linear combinations 

e± = i(4.±4J (120) 

Writing equations ( |117| - |119| ) in terms of these variables, taking the + and — linear combina- 



tions of ( p.l7|) and ( |11^ ), we get 

+ (2r,; + r+)e+ + r_(e-- 03^,3) = (121) 
+ r_e+ + r+(e_- 03^,3) =0 (122) 
a%3 + 2r_e+ + 2r+(e-- 03^,3) - C3r^eg3^(-^) = 0. (123) 
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where 

r±^i(r,^±(r, + r,^)). (m) 

These equations show how the chemical potentials for the other combinations of species, 
namely ^_ and the right handed quark asymmetry are not directly sourced, but develop 
only after the initial generation of a nonzero value for ^-|.. Curiously, the quark asymmetry 
can accidentally vanish if the interaction rate r_ is zero, but this is a special case and not 
generic. We also see that generation of ^+ is damped by the fast reaction rates T^+2Thf. This 
is the suppression which has led other authors to conclude ^_|_ = and to therefore ignore 
the source for this combination. The suppression is numerically not very large, however, 
because the diffusion length of Higgsinos is comparable to the scattering length associated 
with the rates r+ + 2Thf- Thus a significant asymmetry of hi + h2 can build up despite the 
damping term in ( |121|) . 

The Tm term in eq. ( |123| ) makes ( |121| - |123| ) difficult to solve. We can proceed by approxi- 



mating the interaction terms in eq. ( |123| ) to be in thermal equilibrium. Then one can ignore 
the diffusion operator compared to the damping terms, and solve for 

(In fact only the region z > will be relevant for what follows.) The approximation can be 
justified by estimating the relative sizes of Vg^g^ and the rate terms. Below we will show 
that ~ exp(— /c/jx), where ~ T^/6Dh for large r_|_. Thus 

%3 ~ Dgkl^g, ~ r+^e,3 « r+e,3, (126) 

A subtle point here is the use of the large diffusion length /c^^ of the Higgsinos, rather than 
the much shorter one which might be expected for quarks, due to their stronger interactions. 
Physically this occurs because the diffusion tail of comes from the tranformation of 
efficiently-diffusing Higgsinos into quarks, by interactions like h + tn ^ q^. 

Using (|125|) in eqs. ( |121[ - |122|) , the problem is reduced to two coupled equations: 



+ ^(-^) ^- = Sh (127) 

T^Hi^ + J^^ZT^^(-^)^- = -J-TT^^(-^)^+- (128) 
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We note that ^_ is generated only in the broken phase, behind the bubble wall. If the top 
quark helicity flip rate was zero, ^_ would decouple completely and the system would be 
reduced to a single equation for ^_|_, easily solvable by Greens function methods. However 
the helicity flip rate is large, since rrit ~ yVc/Tc ~ 1, where Vc is the Higgs VEV at the 



critical temperature T^. Thus ( 127 - 128 ) cannot be further simplified without additional 
assumptions. One such assumption which is quite plausible is that the ratio 



R ^ T^/T, 



(129) 



is small, in which case ^+ and do approximately decouple. Even with a rather mild 
suppression R < 1/2, we can drop the ^_ and r_-terms in eq. ( |127|) . The solution then 



proceeds in two steps. First one solves eq. (|127|) for which now reads 



This can be done using the Greens function: 

U^) = l r dyG{z-y)SH{y), 



(130) 



(131) 



where 



''+^ X > 



X < 




F = 2F;,/ + F+. 



(132) 



Although ^_ is suppressed relative to by the factor R, one cannot simply neglect it, 
because in eq. ( |125| ) for the contribution from ^_ is enhanced by 1/-R relative to that 
from C,+- In our approximation the equation ( |128|) for ^_ becomes 

r r r r 

-L-i-J-m \^ -L — -L 777. 



2T: + F„ 



2F , + F, 



where 1^+ is the integral solution (|131|) . Defining a scaling parameter 

F. 



7 = 2 + - 



(133) 



(134) 



equation (|133| ) may be written as 



-Rr^e{-z)^, 



(135) 
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which can also be solved by Greens function methods. We only need to know the solution 
in the symmetric phase, because the sphalerons which are being biased by C,q-^ are highly 
suppressed in the broken phase. We find 

e-(^) = -|i?r^ r dyG^{z,y)^4y) , z > 0, (136) 

J — oo 

where the new Greens function is 

G>iz,y) = ^e--^/^4-e(-y)^e-"-^ 

+ e{y)e{z - y) [e"--yl^^^ + +e{y - z) [e--'l^^ + } ^^^'^^ 

(of which only the first term enters in ( p.36| )) with 



It can be seen that the solution ^„ vanishes not only when i?— >0, but also when Fm— >0, as it 
should. Inserting the solution ( |131| ) for i+{y) into ( |136| ) and performing the y integral, one 
can write ^_ as 

Uz) = -l^R^^e-^^-/^-^'-^' / dyg.{y) Sniy), z > 0, (139) 
with the kernel ^_ 

With these solutions for the Higgsino chemical potentials, eq. (|125|) gives that of the 
third generation left-handed quarks. The first and second generation quark potentials are 
determined by eqs. ( |11U| , |112| , [1161 ). We are now ready to consider how the quarks bias 



sphalerons to produce the baryon asymmetry. 
5.1 Baryon asymmetry 

Local baryon production is sourced by the total left-handed quark and lepton asymmetries 
in front of the bubble wall. In the present scenario, there is essentially no lepton asymmetry. 
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Thus the source for baryon production by the passing wall is just the left-handed quark 
asymmetry, ^g^, which enters the baryon violation rate equation as 

^ = ir.p.(f„-^^). (141) 

Here Fsph = KgphCt^T^ is the Chern-Simons number (CSN) diffusion rate across the energy 
barrier which separates A^- vacua of the SU(2) gauge theory, where Kgph = 20 ± 2 [p3| . 
The second term describes sphaleron-induced relaxation of the baryon asymmetry in the 
symmetric phase (more about which below). Using (|110| , |112| , p.l6| and p.25| ) one finds that 



the quark chemical potential created by the classical CP violating force in the wall, combined 
with fast Yukawa and strong sphaleron processes, is 

~ 3 (^gi + + iqs, ) 

= 3C(K,)(i?^+ + ^), ^>0. (142) 



where the factor of 3 counts the quark colors, i+{x) and i-{x) are given by eqs. ( |131| ) and 
( |133|) , and 

with Ki 2 defined analogously to in eq. ( |113|) , and C3 given by (|116|) . The coefficient C{ki) 



encodes the essential information about the effect of the squark spectrum on our results, 
as will be shown below. We remind the reader that ( |142| ) is valid to leading order in an 
expansion in R = r„/r+, which is assumed to be smaller than unity; recall that ^_ is of 
order R in (p9|). 

The second term on the r.h.s. of (|141|) is the Boltzmann term which would lead to 
relaxation of the baryon number if the sphaleron processes had time to equilibrate in front 
of the bubble wall. This would be the case if the bubble wall was moving very slowly. Thus 
the hb appearing here is related to the left-handed quark and lepton asymmetries, fig and 
fii, that would result from equilibrating all flavor- changing interactions which are faster than 
the sphaleron rate in the symmetric phase. Thus A is given by 

A^^fics = 9l^q + Y.f'k^ (144) 

where fics is the chemical potential for Chern-Simons number, and the latter equality follows 
from the fact that each sphaleron creates nine quarks and three leptons. Because of efficient 
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quark mixing, all quarks have the same chemical potential fig. In the leptonic sector however, 
the mixing may be weak (depending on the neutrino and slepton mass matrices) and each 
flavor asymmetry may be separately conserved. To solve for these chemical potentials, one 
must determine which interactions are in equilibrium on the relevant equilibraton time scale, 
which depends on the spectrum of super symmetric particles carrying baryon and lepton 
number. In the usual wash-out computation in the broken phase, the appropriate time scale 
is the inverse Hubble rate, which means that even the feeblest Yukawa interactions leading 
to Ci^-equilibration |3|] are considered to be fast. In this case, using the notation 

na = Na-N-a = na^, (145) 

o 

where Ha = 1(2) when a refers to a fermion (boson), one can show that for the SM 

ns = = ^(6 X 3 X 2)^ ^ = 1 ^ (146) 

and similarly J2if^k = 27^b/T^, which, when inserted in ( |144|) gives the familiar result A = 
13/2. For electroweak baryogenesis, however, the relevant time scale is the inverse sphaleron 
rate in the symmetric phase, and therefore none of the right-handed leptons will have time 
to equilibrate (r is in fact a border-line case with chirality flipping rate comparable to the 
sphaleron rate; but we take it also to be out of equilibrium). Then, with the SM spectrum, 
which would apply if all squarks were heavy, J2i f^k = 3nB/T^ and hence A = 15/2. If there 
are Ngg flavours squarks which are light enough to be present at T = 100 GeV, one has 

A = l(l + ^) ' + 3. (147) 



2 V 6 

It is straightforward to generalize A to the case of an arbitrary number of light left-handed 
sleptons, but the expression is cumbersome because of the multitude of possible mixing 
scenarios in the leptonic sector, and we omit it here for the sake of simplicity. 

Moving to the wall frame, the time derivative in ( |141| ) becomes — v^jdz, and it is 
easy to integrate the equation to obtain the baryon asymmetry: 

3 Tsph 

where 



/ dzU{z)e~'-\ (148) 
Jo 

■iA Fsph 
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(149) 



The integral over z in (|148|) can be done analytically. The baryon-to-entropy ratio, 



?7a = n^jn^ ~ Iub/s, can then be written as a single integral over the source function 
SH{y)- 

= 945^EPh4l C{K.) R r dy [g^{y) - ^ , Sniv), (150) 

Sn^v^g^ J-oo \ Q-fa^iv^ + DhkB) J 

where we have scaled the variable y to units 1/T, and the new kernel, arising from performing 
the 2;-integral over ^+{z), is given by 

n-i f / pk-y / 1 1 \ \ p'^+y ) 



A;+ — A;_ 1 \ A;_ + /c^ \k+ + k^ + ks / / k+ + k 



B 



(151) 

The ratio of the contributions coming directly from ^+ and from the indirectly sourced ^_ 
is controlled by the parameter Tm/ {a^^v^ + Dfiks)). Using typical values for the other 
parameters, the ^_ term turns out to be significant for ^ 0.01, and subdominant for 
larger or smaller wall velocities. 

5.2 Sources 

We still need to calculate the source 5*^ appearing in the above equations. It is given by the 
thermal average ( pHD where the force 6F corresponds to the CP-violating part of the classical 
force, eq. (0). For a Higgsino of helicity A, using kj^ = 1 since Higgsinos are fermions, we 
have 

A /\Pz\ \ I 2ai\ii nc;o^ 

where the sign ± is defined to be the sign of — |m2|, since the lighter (heavier) of the 
local mass eigenstates is the Higgsino-like particle when /i < m2 (/i > vn-i). The absolute 
value on \pz\ comes from the relation between spin and helicity: sp^ = \\'Pz\- The average 
(fpj is very accurately approximated by the fit 

^ X%. + 6X± + 2 

where x± = m±/T, and using Maxwell-Boltzmann statistics, one can show that 

/M\ = a - x±)e-^^ + xlE^{x^) 
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where K2{x) is the modified Bessel function of the second kind and Ei{x) is the error function 



m 



In deriving ( |152|) , we imphcitly assumed that the charginos are hght compared to the 



temperature. In the limit that they become heavy they must decouple however, and as 
a result the damping rates r± for Higgsinos to be transformed into quarks/ squarks, in eq. 
( |123|) , must go to zero. The approximations (|125| ) and ( |126|) would consequently break down. 
In an exact treatment, one should solve the equations ( p.21| - |123D numerically in these cases. 



We will instead adopt a simpler approximation, incorporating the effect of decoupling on 
the chargino source Sh with a suppression factor n(m^)/n(0), that is, the ratio of thermal 
densities for a particle of mass relative to a massless particle. In this approximation, Sh 
becomes 

SH,eS = - ^±)e-'^ + 4^1 K^i)"- (155) 

The effect of this modification is small for chargino masses up to 200 GeV. Beyond this it be- 
comes crucial for suppressing baryon production from particles too heavy to be present in the 
thermal bath. For ^ 200 GeV, our results should be understood to have a multiplicative 
uncertainty of order unity arising from this approximation. 

To fully specify the source term, we must also give the functional form for the spatial 
variation of the Higgs field condensate, since this enters the mass eigenstates 'm± and the 
CP-violating phase 6'^ in the above formulas. Because our source is proportional to the 
combination HiH'2 + H2H[ of the two Higgs fields, our results are not sensitive to changes 
in the ratio tan/? = H2/H1, which is in fact known to be nearly constant for bubble walls in 
the MSSM — at least for generic parameter values, including those that give a strong enough 
phase transition. This is in marked contrast to other analyses where the source was assumed 
to be proportional to H1H2 — H2H[. It thus suffices for us to use a simple kink profile 

Vr 1 



u{z) ^ g^Hl + Hl = g-^-{l- tanh ( - ) ) (156) 



with ui = Msin/? and U2 = ucos/3. Here £^ ~ 6 — 14/T [ Pq , |T2[ is the wall width and Vc 
is the value of the Higgs condensate at the critical temperature. This VEV has the usual 
normalization in the vacuum vt=o — 246 GeV, while the requirement for a strongly enough 
first order transition (to avoid washout) is Vc/Tc > 1.1. 

We display the source ( |155|) as a function of position relative to the wall (at z = 0) in 
Fig. 1(a), using the parameters values 1712 = 150 GeV /i = 100 GeV, 5^ = arg(m2;u) = tt/2 
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Figure 1: (a) The source for baryogenesis from the chiral classical force, eq. ( |155| ), for the 
parameters /i = 100 GeV and m2 = 150 GeV and = 10/T (solid line) and = 14/T 
(dashed line), (b) The left-handed quark asymmetry ^g^, eq. ( |142| ), for the same parameters. 



The distance from the center of the wall z, is measured in units 1/T. 

and Vyj = 0.3 for two different wall widths: i^j = 10/T and = 14/T. For the parameters 
left unspecified above we use the following standard reference values: 



20/T Thf = 0.013 r„ = 0.007T r+ = 0.02T, R = — = 0.25 



r. 



10/T, t;, = 120, T = 90, tan/5 = 3. 



(157) 



In the limit that u"^ ^ /i^, the source would be a symmetric function of z since it 
would be proportional to {u^)'" . However for finite /i and m2, the actual 2;-dependent mass 
eigenvalues appearing in the coefficient of (u^)'" depend on rather than its derivatives, 
hence the departure from the symmetric form. In figure 1(b) we plot the profile for left- 
handed quark number, for the same set of parameters. As expected, the spatial extent 
of the quark asymmetry is roughly the diffusion length of the Higgsinos, Dh/vw ~ 60/T (see 
the remarks below eq. (|126| )). For a smaller wall velocity the diffusion tail extends further, 
but the amplitude of in the tail gets smaller because then the damping has more time to 
suppress the chargino asymmetry. 
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The rates quoted in (|157|) are rough estimates, obtained from an approximate computa- 
tion of a subset of relevant 2^2 reaction rates, and higgsino decay rates, when kinematically 
allowed. For example h2W^^ijiqL, gives 



aw 



fe + >n-r). (158) 



where s ~ 20T^ is the center of mass energy squared and m?^ ~ mfj + 0.9T^ is the thermal 
mass of the right handed squark . The soft SUSY breaking mass parameter mfj is taken 
to be negative, mfj ~ —60^ GeV^, as indicated by the need to get a strong enough first 
order phase transition |ll|. The rates (|157|) correspond to a conservative overestimate by a 
factor of 5 over the total averaged contribution from various scattering channels (for how to 
perform the thermal averages, see ref. [^). The decay rates have a fairly strong dependence 
on higgsino mass ~ /i. For example /i2L~^^i?&i gives 



^a^Xni,^,^), (159) 

1 e_ V ' mi ' mi ^ ' 



where A(x, z) = [x — y — z)^ —4:yz and m^^ ~ 0.76T |^. For small /i ^ 130 GeV, the decay 
channels are not open, whereas for large enough /i they dominate over the scattering con- 
tribution. However, our numerical results for rj^ are fairly insensitive to changes in various 
rates; for example decreasing F^ by a factor of 5, appropriate for /i <, 130 GeV would in- 
crease riB by about 30 per cent, whereas incrasing it by a factor of 5, appropriate for fi ~ 500 
GeV, would decrease r/s by about 40 per cent. This scaling is somewhat weaker than the 
naively expected rjB ~ l/v^TV dependence, because the damping effect due to faster rates 
is initially being compensated by more efficient transport from the chargino sector. (We 
have checked that the naive scaling eventually follows for values of F+ large enough that the 
transport effect has been saturated.) Nevertheless, the relative insensitivity of the results on 
the rates warrants our use of the rough estimates (p. 5 7]) in our numerical work. 



6 Results 

Dependence on squark spectrum. Let us first consider the dependence of rjg on the 
squark spectrum. This is contained in the parameter C{ki), some representative values of 
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I'll 7 

light squarks 




A 11 

All 


u 


All R-chiral 





All 3rd family 





tL, II and 


2/41 


II and 


3/16 


Ir and hp^ 


3/8 


tR only 


10/23 



Table 1: Multiplicative factor C(/«j) containing the dependence of rjg on the squark masses 
for particular choices of the light squark spectrum. 



which are given in Table 1. For certain choices of squark masses, C{ki) = 0, which re- 
flects the approximation we made of taking the strong sphalerons to be in equilibrium; it 
is well known that these interactions tend to damp the baryon asymmetry if, for example, 
no squarks are present |3^. In these cases the baryon asymmetry is not really zero, but 
comes from l/Tgs corrections which we have not computed. Ignoring such corrections, one 
sees the clear preference for the minimal possible number of light squark species from C{ki). 
This is fortuitous because it coincides with the need for a single, light, right-handed stop 
in order to get a strong phase transition. If the left-handed stops and sbottoms are also 
light (which, incidentally, is incompatible with the large radiative corrections needed for the 
Higgs mass to satisfy the experimental lower limit, as well as rho parameter constraints) the 
baryon asymmetry is reduced by a factor of ten. Thus, considerations both of the initial 
baryon production and the preservation from washout favor the "light stop scenario." Since 
the effects of the spectrum are trivial to account for in the final results, being just an over- 
all multiplicative factor, we shall henceforth concentrate only on the most favorable scenario. 



Velocity dependence. The dynamics of the phase transition, even apart from CP- violating 
effects studied here, is a very complicated phenomenon, involving hydrodynamics of the fluid 
interacting with the expanding walls, and reheating effects due to the latent heat released in 
the transition Although the originally spherical bubbles quickly grow and reach some 
terminal velocity, inhomogeneities can subsequently develop. This occurs when the shock 
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waves from the bubble expansion heat the ambient plasma and thereby reduce the latent 
heat released as regions of space are converted from the symmetric to the broken phase. 
There is a subsequent decrease of pressure driving the expansion, and depending on model 
parameters, may lead to significant slowing down of the walls. The process of heating by a 
collection of shock waves causes local variations in the temperature as well as fluid velocities, 
with consequent deformation of the shape and speed of the wall. These variations occur on 
the macroscopic length scale of the bubble radius, which is many orders of magnitude greater 
than the microphysical scales that have been discussed here so far. In this sense, eq. ( |148| ) 
gives only the local baryon number at a given position in space after the wall passes by. The 
presently observed asymmetry should be computed by averaging over a region which is large 
compared to the bubble size at the time the phase transition completes: 



where 7]^ is considered as a functional of the locally varying wall velocity. Only if the phase 
transition is very strong, so that there is a high degree of supercooling, will the reheating 
effects leading to inhomogeneities be small or negligible. 

In addition to the possibility that Vw has spatial inhomogeneities, it is also interesting to 
study the dependence on simply because its value is not yet known with great certainty, 
although some progress has recently been made |2^, Our treatment takes into account 
the back-reaction effect on the baryoproduction (washout by sphalerons), so our results are 
valid for arbitrarily small wall velocities. In Fig. 2(a) we plot tjiq = rjB x 10^^° as a function 
of for /i = 100 GeV and ms = 50, 100, 150 and 200 GeV, and in Fig. 2(b) for four different 
values of the wall width i^, with fi = 100 GeV and m2 = 150 GeV. 

The peak occurring at ~ 0.01 for some parameters in Fig. 2 (a), first observed in 
|20| , is due to the contribution from the Q- term in ( |150| ). This is enhanced by a factor 
(fu, + DhkB)'^, which for the assumed parameter values peaks near Vw ^ y/Dfiks — 0.01. 
Because of the back-reaction, the baryon asymmetry vanishes when the wall velocity goes 
to zero. The peak is prominent only for the values of m2 ~ /i however, and the typical 
velocity depencence of rjio is not quantitatively very large as a function of velocity. It is quite 
complicated however, in that for special parameter values the asymmetry can accidentally be 
small or zero. The crossings through zero arise as follows: for relatively large Vw the baryon 
production in the diffusion tail dominates over the opposing contribution generated near the 
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Figure 2: rjiQ as a function of wall velocity for yU = 100 GeV and sin 5^ = 1 for: (a) a varying 
gaugino mass parameter m2 = 50, 100, 150 and 200 GeV and: (b) for a varying wall width, 
C = 6/T, 10/T, 14/T and 18/T. 



wall (see the generic form of the distribution in Fig. 1 (b)). For small wall velocities 
the length of the diffusion tail increases as D/vw, but the amplitude of the asymmetry gets 
smaller due to interactions, which have more time to damp the asymmetry. Moreover, the 
contribution from the part of the diffusion tail extending beyond l/ks is cut out, because 
the baryon asymmetry is already relaxing due to sphaleron washout beyond that distance. 
As a result the contribution from the tail eventually becomes the smaller one, leading to a 
cancellation between the two contributions that give the net asymmetry. 

While the uncertainty in Vy^ at present is not necessarily the dominant one for estimating 
the baryon asymmetry, determining t]iq to high precision for a given set of chargino mass pa- 
rameters would need careful hydrodynamical modelling of the bubble wall expansion. Also, 
even rather small fluctuations in rjB can have interesting consequences elsewhere: for ex- 
ample they can seed the generation of large fluctuations in leptonic asymmetries in certain 



neutrino-oscillation models pO| with potentially large effects on nucleosynthesis. 



Dependence on chargino mass parameters. The most important supersymmetric in- 
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Figure 3: Contours of constant CP- violating phase 5^, corresponding to baryon asymmetry 
?7s = 3 X 10~^° for (a) = 0.1 and (b) = 0.01. Mass units are GeV. 



puts directly affecting the baryon asymmetry are the chargino mass parameters m2 and 
and the CP-violating phase 6^ = arg(m2/i). In Fig. 4a and 4b we plot the contours of con- 
stant \6fj_\ giving the desired baryon asymmetry 77^ = 3 x 10~^° |^l|] in the (m2, /i) plane. The 



baryoproduction is most efficient for small masses, m2,/^ ~ 100 GeV, but baryons can still 
be copiously produced for \m2\ and ~ 500 GeV and higher. In the best large 
enough baryon asymmetry can be produced even with a very small explicit CP-violating 
angle of order a few xlO^'^, comfortably within the constraints coming from electric dipole 
moment searches [Q. 

7 Conclusions and outlook 

We have presented a detailed analysis of electroweak baryogenesis in the minimal supersym- 
metric standard model (MSSM) using the classical force mechanism (CFM). We argued that 
the dominant baryogenesis source in the MSSM arises from the chargino sector. We also 
commented on a recent controversy regarding the parametric form of the source appearing 
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in the diffusion equations. The resolution is that all different formalisms agree with the 
parametric form; however previous authors neglected the particular source considered here 
for the linear combination of Higgsinos Hi + H2, on the grounds that it is suppressed by 
the top (s)quark Yukawa interactions. We have shown that this suppression is quite modest, 
a factor of order unity, which is much milder than the intrinsic suppression suffered by the 
competing source Hi — H2, due to the near constancy of the ratio H1/H2 throughout the 



bubble wall p6|, p 



Our present work differs in several ways from our earlier published results using the CFM. 
First we have presented a treatment in terms of the physical, kinetic variables characterizing 
the WKB states rather than in terms of canonical variables, which are gauge dependent, and 
in terms of which the recovery of a gauge independent physical result is not always trans- 
parent. While this is mainly a matter of (considerable!) convenience, there is also a slight 
physical difference in our results due to a slightly different form for equilibrium ansatz corre- 
sponding to each parametrization (see appendix B). We noted that a definitive determination 
of the correct form will have to await the outcome of a more fundamental computation, as 
will the correct treatment of 'spontaneous' baryogenesis, which rely completely on the rele- 
vant form of the collision integral. Much more importantly to quantitative changes in our 
results is that in our treatment in |jl8[ we misidentified the sign of the hypercharge of one of 
the Higgsino states. This prevented us from realizing that the top Yukawa interactions tend 
to damp the appropriate combination of Higgsino currents, Hi + H2- Here we developed a 
new set of diffusion equations where this effect is treated correctly. Thirdly, here we have 
considered several different choices for which flavors of squarks are light compared to the 
temperature, and found that the one adopted in |jl8|] (all squarks light) is among the less 
favoured possiblities for baryogenesis, because of strong sphaleron suppression. 

We have given a complete derivation of the CFM formalism, starting from the basic 
assumption that the plasma is adequately described by a collection of WKB-quasiparticle 
states in the vicinity of a varying background Higgs field. We derived the dispersion relations 
for squarks and charginos, and showed how to identify the appropriate kinetic momentum 



variable, the physical group velocity, and the force (see also |^). We pointed out that the 
force term for the current combination Hi — H2, obtained in a recent publication [^], is 
absent in physical variables and hence this current is not sourced in the CFM. This is a 
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very sensible result, because the term in question can always be removed by a canonical 
transformation, or equivalently, a field redefinition. We also derived the diffusion equations 
and the source terms appearing in them starting from the semiclassical Boltzmann equations 
for the quasiparticle states. 

We have studied the baryon production efficiency in the MSSM as a function of the pa- 
rameters in the chargino mass matrix, the wall width and the wall velocity. The dependence 
on wall velocity is rather complicated, and intertwined with the dependence on the chargino 
mass parameters; the generated asymmetry generically changes sign as a function of f^, but 
the value of v^, where this crossing takes place, and the the functional form of risivw), are 
quite dependent on mass parameters. However, for large regions of parameters an asymmetry 
rjiQ = lO^^riB of the order of several hundred could be created, implying that a CP-violating 
angle of arg(m2/i) ~ 10~^, and in best cases even arg(m2yu) — a fewxlO^'^, suffices for pro- 
ducing the observed asymmetry of r]io ~ 3 [^. Such small phases are consistent with the 
present limits from the neutron and electron dipole moment constraints . 

We finally emphasize that our formalism disagrees in detail with various other meth- 
ods of computing the source in the diffusion equations. This is particularly significant with 
regard to references |^4[ 0, |2^, [1^, which all claim to be valid in the thick wall regime, 
where our method was designed to work. This is troubling, because one expects that differ- 
ent methods should agree when the same physical limits are taken. In particular, we have 
shown that classical force mechanism does not give rise to a source of the parametric form 
~ H[H2 — H1H2, found by references [^, |2^, |17|, |22[ for both squarks and charginos. We 
do not know a definite solution to this problem, but a possible origin for the discrepancy 
could be that the methods p^ , [1^ |2^, perform an expansion in the mass, or the vac- 
uum expectation value, divided by temperature (the mass insertion expansion [^]), before 
taking the gradient expansion. In the WKB approach on the other hand the background 
is treated in a mean field approximation and one performs the gradient expansion around 
this classical background. In other words, in the WKB-picture the mass insertion expansion 
has been resummed to infinite order before the gradient limit is considered. While one can 
formally expand the CFM-source resulting from a WKB-analysis in mass over temperature, 
one should in general not expect that taking these two limits is commutative. In partic- 
ular, quantum refiection is completely absent in the WKB approximation, but is certainly 
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present in the mass insertion expansion. The issue of how to properly account for both the 
semiclassical and the quantum effects, or to interpolate between them, is certainly worth 
further study, and some published results from a work aiming to a derivation of appropri- 



ate semiclassical Boltzmann equations from first principles can be found in references |32, 31 
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Appendix A: Collision Terms in Linear Expansion 

We show here how the collision integral on the r.h.s. of the Boltzmann equation gives rise 
to the terms damping the perturbations from equilibrium. For illustration, let us consider a 
two body process with ingoing WKB states i with four momenta Pi (with pn corresponding 
to the distribution on the l.h.s. of the Boltzmann equation) and outgoing WKB states / 
with four momenta pf-. 

c[f,] = ^ / \M\\2nrd'(j:pi)nf)] (lei) 

Jpi2,Pf I 

where is the matrix element calculated between the WKB states, to first order in 

derivatives of the background, Jp means / d'^p/2£'(27r)'^ and the statistical factor V[fj] is 
given by 

Vm = UUl T //i)(l T //2) - //i//2(l T /.i)(l T h2). (162) 

We only attempt to compute the collision integrals to the zeroth order accuracy in gradients, 
so that the integral measures and 5-functions are the same as in the usual fiat space-time 
considerations. Most part of the derivation consists of manipulating the statistical factor 
Inserting the ansatz ( [78D to ( |162| ) and expanding to the first order in 5f one gets 



\ In Ji2 Jfl Jf2 I 
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where /j(/ii)'s are the distributions given by the first term in the ansatz (|78D and /^'s are 
the unperturbed distribution functions. V[fi{fii)] is nonzero only for inelastic scatterings, 
whereas all reaction channels create nonvanishing collision terms proportional to 6fj; let us 
consider these terms first. 

The entire collision integral corresponding to 6f j-teims in ( |163|) can be written as 

C[6fi] ~ 6Upi)r{pa) + fiPii) f 5/.2fe2)a(pa,P,,2) 

-E/ KiPfn)Gfn{p^l 5 Pfn) 
n •'Pfn 

= Sfa{Pi)'rM-SF{pa) (164) 

The abbreviated notation here is used to highlight the fact that only the first term is di- 
rectly proportional to Sfn^pn), whereas all the others contain smeared integrals over 6fj- 
distributions; the exact forms of the functions a{p,k) and Gi{p,k) are not relevant for us. 
T{pii) on the other hand is the usual thermally averaged interaction rate, which, neglecting 
the Pauli-blocking factors, is given by 

mi) = I ^,fUP^2KVr.la,^J). (165) 

where 



Pf I 



is the invariant cross section for the process i-^f multiplied by the invariant flux (see for 
example p9|). 



The 6f F-term in (|164| ) clearly causes damping away of kinetic fluctuations, with a mo- 



mentum dependent relaxation scale given by the inverse of the rate ( |165| ) . The "noise term" 
6F physically represents the process of further thermal redistribution of the states which goes 
on alongside the relaxation of 6fi to zero. Thes are random processes which occasionally op- 
pose the relaxation process. However, while the integrated over pn moments of 6F{pii) are 
comparable to the moments of the first term, their naive inclusion to the moment equations 
would be incorrect, since kinetic relaxation depends sensitively on the shape of the entire 
distribution function. For example, the condition {{pn/ EiQ){6fT — SF)) = would lead to 
vanishing of the (kinetic) relaxation term for the velocity perturbation in moment equa- 
tions, whereas in reality the kinetic relaxation process is halted only if the collision term 
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( |164|) vanishes identically for all momenta. The effect of the noise terms is further reduced 
by the fact that in all elastic channels adding the scatterings from particles and antiparticles 
tend to cancel the noise part, while the contributions to the relaxation terms are equal and 
add. In the diffusion approximation then, the first moments of the part of the collision term 
containing (5/j's are given by 

{C[Sf,]) ^ 

{v.MSfj]) ^ {vpMiiPz))^, (167) 

where the average (■) is as defined in equation (|87D . In the case of first moment we have also 
assumed that T{p) has only a weak momentum dependence so that it can be replaced by its 
thermal average T. This is of course the place where we implicitly truncate our momentum 
expansion to the first two terms, and it should be a very good approximation. Adding up 
the contributions from all possible channels, one obtains the result (0). 

Let us next consider the V[f{fij)]-pa.Tt of the statistical factor ( |163|) . Expanding to the 
first order in fi/s one finds 

nfi^)] ^ -M2{u + ^^2 - Oi - 02), (168) 

where ^ = fi/T and we also neglected the Pauli blocking factors in the final states. This 
expression obviously vanishes for the elastic channels, whereas for inelastic channels it gives 
the contribution 

c[^^,]^f^MUp^l)J:^J (169) 

j 

where Ti^pn) is an expression analogous to ( |165| ) and is the signed sum over the 

chemical potentials appearing in ( |168|) such that the term has a positive sign. The first 
moments of the inelastic collision term in ( |169| ) then are 

{v,Mf^,]) ^ (170) 

where 

f. = / f^^Up,,)/Na (171) 

with Nil = I d^p/j'o (see equation (^)). Adding up all inelastic channels affecting a given 
species i, one arrives to the equation (BTI). 
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Appendix B: Equilibrium ansatz in canonical variables 



Gauge invariance. In previous treatments 0, 18| an ansatz different from (|77|) was adapted 
for tlie local equilibrium function: 

f(pc, x) = —T, — ^- — M , , + Sf(p, x), (172) 

where Pc is the canonical momentum. A technical problem with the canonical variables is that 
both Pc and fi are phase reparametrization, or "gauge" variant quantities. This can be seen 
by observing that any physical quantity (e.g. local number density) is obtained integrating 
over the momenta. The integration measure d^pc is unchanged by gauge transformation 
Pc ^ Pc + o^cp, so that a system with fixed number density is described by a different value 
of /i in two different gauges. 

In the previous treatments ^ and equations for physically meaningful quantities 
were recovered using the condition that the system be unperturbed from equilibrium far in 
front of the wall (at z oo) i.e. 

/ d^Poof = [ d^p^—T, — ;— r i Mnhvs ^ as 2:^ 00 (173) 

where poo is the physical momentum at infinity. For the case of a fermion with complex mass 
discussed in section 2, Eqns. (0) and (|18D give 

s9' 

Poo = Pc,oo + Sep— - acp (174) 

(since m — > 0), where acp = a' + ScpO' /2 in the leftchiral sector for example. One then 
identifies the physical chemical potential as 

s9' 

Aiphys = + Vw'^w{scp^ acp). (175) 

The equations are then most conveniently rewritten in terms of /iphys and poo, and solved 
with with the boundary condition /iphys = at +oo. Note that in both and fl^ a specific 
gauge was chosen, which can be read off from the dispersion relations adapted as acp = 
in P], and a^p = Scp^ in |jl8[. While in the former case a transformation from the original 



canonical variables had to be performed (cf. section 4, page 2962 in 0), the implicit gauge 
choice of |jl8| required no such transformation, since /iphys = /i in this gauge. Note in par- 
ticular that the terms proportional to the linear combination of scalar fields H[H2 — H2H1 
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appearing in canonical momenta (1^), (|6^ and (|69|) , can entirely be absorbed into the redef- 
inition of /iphys, so they will not provide new sources even when treating the problem using 
the canonical momentum. 



Comparision of the ansatze. Making use of ([20| ) and the dispersion relation (|T9D one can 
show that the canonical and kinetic momenta are, to linear order in 9' , related by 

s9' 

p^{l-Scp—)=Pc-Oicp, (176) 



where cj = ^Ju"^ ~ ■'■^ kinetic variables the ansatz ( |173| ) may then be written as 

= g/3b.(..+..pUp..s] ± 1 + ^f^P^ ^) + /?^-7«.ScP^^/', (177) 

where /' = (1/e^ ± 1)' {x = (3uj). One can thus identify the physical chemical potential yUphys 
in the canonical variables with the chemical potential for our physical WKB-quasiparticles 
appearing in the ansatz (^). The chemical potentials in ([78D and ( |172| ) are thus only 
separated by an unphysical gauge-transform. The distributions differ however, by a term 
that cannot be transformed away: to make ( |1 72| ) completely agree with (^), one should 
have 

5f = 6f- f3v^j^sc/-^'^f. (178) 

The latter term does not vanish in equilibrium however, so the two ansatze do correspond to 
two physically different equilibrium conditions. The difference is only nonzero in the region 
of the wall however, as it vanishes when |m|— ^0. Including this term in our Boltzmann 



equations would contribute to source term, making it equal to the one used in [ffH . 

It is clear that the difference between ansatze corresponds to which energy momentum 
- canonical or kinetic - is the appropriate one to take as that conserved in the local inter- 
actions between particle states modelled there. We have argued in the main text that the 
kinetic momentum has the more direct physical interpretation, and this argument is backed 



up by results from a more sophisticated treatment [|3ll , |32| . However, a complete derivation 



of the transport equations using the formalism of ||3T|, |3^ will be needed to settle the issue 
unambiguously, while in practice the numerical results are not particularily sensitive to the 
difference. 
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